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Abstract 

In this paper, a distributed output regulation problem is formulated for a class of uncertain nonlinear multi-agent systems 
subject to local disturbances. The formulation is given to study a leader-following problem when the leader contains unknown 
inputs and its dynamics is different from those of the followers. Based on the conventional output regulation assumptions and 
graph theory, distributed feedback controllers are constructed to make the agents globally or semi-globally follow the uncertain 
leader even when the bound of the leader’s inputs is unknown to the followers. 
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1 Introduction 

The past decade has witnessed a rapid development 
in the field of multi-agent systems and fruitful results 
have been achieved for the leader-following coordination 
problem. In recent years, distributed output regulation 
of multi-agent systems has been proposed to provide 
a general framework for leader-following consensus in 
linear or nonlin e ar cas es (I Hong. Wang fc Jiang 1. IM1 
iDong fc Huand . l2Q14ll . Internal model approach, de- 
veloped to solve the conventiona l outpu t regu lation 
(| IsiYri. Marconi fc SerraiiH . 120031 : iHuand . 1200^ , was 
effectively used for the distributed design, especially 
for nonlinear agents. For example, a cooperative out¬ 
put regulation problem was considered for a class of 
nonlinear uncertai n multi-agent syste ms with unity 
relative degree in ISu fc Huand (|2013ll . while agents 
in the output-feedback form were co nsidered with an 
autonomous leader and no-loop graph (|Ding I . l20I3ll . 
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It may be restrictive or unpractical if we always consider 
an autonomous leader without unknown inputs, espe¬ 
cially in the case when the leader is an uncooperative tar¬ 
get or contains unmodeled uncertainties. Therefore, it is 
necessary to study multi-agent control when the leader 
contains (unknown) inputs. In fact, the generalized out¬ 
put regulation (GOR) problem to track an exosystem 
(or a leader) with extern al inputs was discussed in many 

S ations including | S aber i. Stoorvogel &: S annutil 
and I Ramos. Celikovskv fc Kuceral (l2004l) . On 
the other hand, a distributed problem was investigated 
when the agent dynamics are double integrators to track 
a leader with an un known but bounded acceleration 
in ICao fc (1201 21). and then a similar design was 
given in Li. Liu. Ren fc (|2013D by assuming that 
the leader and followers share the same dynamics. To 
our best of knowledge, there are no general results on 
nonlinear multi-agent control when the unknown-input 
leader and the followers have different dynamics with 
external disturbances. 

The objective of our paper is to study distributed out¬ 
put regulation for leader-following multi-agent systems 
with an unknown-input leader, whose dynamics is non¬ 
linear and may differ from those of the followers. The 
contribution of the work is at least twofold: 

• We extend the distributed output regulation to the 
case when the leader contains unknown inputs and has 
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a dynamics different from those of the non-identical 
followers with (unbounded) local disturbances, and 
provide distributed controls to solve this problem in 
different cases. The results are consistent with the 
existing output regulation results when the leader 
does not have un known inputs and the disturbances 
are b ounded fe.g. lDong fc Huaiia . l20l4lSu &: Huand . 

[Ml). 

• We extend the conventional GOR to its distributed 
version for multi-agent systems with an unknown- 
input leader. Moreover, both global and semi-global 
results are obtained for nonlinear agents with unity 
relative degree, while only local results were obtained 
for a c onventional (single-agent) case in iRamos et al. I 

lIMil . 

Notations: Let R" be the n-dimensional Euclidian space 
and = {s S R” | —M < Si < M, i = 1, ..., n} for a 
constant M > 0. For a vector a;, ||a:|| (or ||a:||oo) denotes 
its Euclidian norm (or infinite norm). diag{6i,.. 
denotes an n x n diagonal matrix with diagonal elements 
(i = 1,.. .,n); col(ai,.. .,a„) = [af,.. for col¬ 

umn vectors Ui {i = l,...,n). A continuous function 
a: [0, a) [0, oo) belongs to class /C if it is strictly in¬ 
creasing and a(0) = 0; It belongs to class /Coo if it be¬ 
longs to class /C with a = oo and lims_>oo Q:(s) oo. 


considered to solve an output co nsensus probl em for the 
exosystem without any inputs in iDing I (|2013il . 

The interaction topology among these agents can be de¬ 
scribed by a graph Q = (V, 8), where V = {0,1,..., n} 
is the set of nodes and £ is the set of arcs. {i,j) de- 
notes an arc leaving f rom node i and entering node j 
(iGodsil fc Rovlell2001ll . A walk in graph Q is an alternat¬ 
ing sequence iieii 2 e 2 - • -Ck-iik of nodes ii and arcs = 
(*m, im+i) G f for / = 1, 2,..., k. If there exists a walk 
from node i to node j then node i is said to be reachable 
from j. Define the neighbor set of agent i as Afi = {j : 
(j,i) S 8} for z = 1,..., n. A weighted adjacency ma¬ 
trix of G is denoted by A = [oy] € R("+i)^fo+i), where 
an = 0 and Qij > 0(ay > 0 if (j, z) € 8). A graph is 
said to be undirected if a^- = aji (i,j = 0,1,..., n). The 
Laplacian L = [kj] £ IR("+i)x(”+i) of graph Q is defined 
as la = hj = —o,ijij ^ *)■ Denote Q as 

the induced subgraph of Q associated with the n follow¬ 
ers. The following assumption has b een widely used in 
coordination of multi-age nt systems (j Hong. Hu fc Gad . 

IM^lSu fc 

Assumption 1 The leader (node 0) is reachable from 
any other node ofQ and the induced subgraph Q is undi¬ 
rected. 


2 Problem Formulation 

Gonsider a group of n -I- 1 agents with one leader (re¬ 
garded as node 0) as follows: 

v=p{v)-i-q{v)w{t), yo = r{v,p) (1) 

where v £ R"" is the leader’s state, and w{t) £ R"“, go € 
R are its input and output, respectively. Here p £ R"'* 
is an uncertain parameter vector, and w(t) is continuous 
satisfying ||zc(/)||oo < I with a constant I > 0. The other 
n (non-identical) agents are followers described by 

(zi = fi{zi,yi,p) ^2) 

lyz = gi{zi, Vi, p) + biUi -\-di, i = l,...,n 

where Zi £ R"*^, z/i € R, di G R, bi > 0. Without loss 
of generality, we take bi = 1 and assume all functions fi, 
Pi, p, q, r are smooth with /i(0, 0 , p) = 0, 5 i( 0 , 0 , p) — 0, 
p{0) = 0, r(0, p) — 0. di is the local disturbance of agent 
z governed by 

uJi = SiUJi, di = Di{p)uji. ( 3 ) 

As usual, we assume Si £ h as no eigenvalues 

with negative real parts (|Huang| . [200411 . 


Given a communication graph Q, denote H £ R"^" as 
the submatrix of the Laplacian L by deleting its first row 
and first column. By Lemma 3 in iHong et al. I (|2006ll . H 
is positive definite under Assumption 1. The distributed 
control law can be constructed as follows: 

Ui = ki{£^i,yi z/j,j £ Afi)i 

— hi (^ 2 , Pi Pj , j £ Afi ) 

where (i £ R"5i with a nonnegative integer n^. and func¬ 
tions ki(-), hi{-) to be designed later. 

To handle this nonlinear multi-agent system with an 
unknown-input leader, we formulate the problem as the 
distributed generalized output regulation problem or sim¬ 
ply distributed regulation problem. It is said to be (glob¬ 
ally) solved for systems (I)-(3) with a given graph G, 
if we can find a distributed control law (4), such that, 
for any (z,(0), z/,(0)) e R"-z+i, ^ g R"., ^,(0) g R-«s 
z;(0) g R”*’,a;i(0) £ R^^z, the trajectory of the closed- 
loop system, composed of (I), (2), and (4), is well-defined 
for all / > 0, and moreover, 

limei(/)=0, ei=yi-po, z = I,...,n. (5) 

t—^OO 

Our problem is said to be semi-globally solved if, for 
any given M > 0, we can find a control law (4) with 
a constant M > 0, such that, for any initial condition 
(z.(0),y.(0)) g R^-+\ p£Mf;jf, ^.(0) g R^^S z;(0) g 


Glearl y, the first-order nonlinear agent in ! Liu. Xie. Ren^ 
(1201311 is a special case of (2), and system (2) was also 


f and Wi(0) g R^b the trajectory of the closed-loop 
system is well-defined for alH > 0 and (5) holds. 


2 



































Remark 1 When n = 1, o ur p roblem becomes GO R 
studied ir ASaberi et al. I ^200 A) and \Ramos et al. I Il200 A) . 
Here we seek non-local output feedback control for non¬ 
linear systems o f the form (2) , while only local results 
were obtained in \Ramos et al. I li200A ) requiring the ex¬ 
osystem’s state. 

Remark 2 Because of the unkno wn expression or typ e 
of w{t), adaptive IM discussed in \Su & Huand \20lA} 
fails t o solve our problem ev en with a time-varying one 
like m l Yana & HuanA \201(\) . In fact, our problem can he 
viewed as a distributed version of GOR to handle those 
exosystems (or leaders) with unknown inputs, which cer¬ 
tainly extends the existing multi-agent output regulation 
formulation when t he leaders have no unknown inputs 
fHona et al. I . \20Hi : liSit & Huand. \20Hi} . 

The followin K assumption was use d for GOR of nonlinear 
systems (see lRamos et al. 1 . 1200411 . 

Assumption 2 There exist two class 1C functions ao(') 
and 7o(-) such that 

IKOII<«o(IKo)||) + 7o(IKt)l|oo). 

Clearly, v{t) is bounded by bi = ao(||r'(0)||) + 7o(0 and 
&2 = max^gRr^K {^g(?;)} is well-defined under this as¬ 
sumption. Moreover, system (1) is Lyapunov stable at 
= 0 when w = 0 with the neutrally stable exosystem 
for nonlinear output regulation as one of its special cases. 

Similar to the output regulation problem, the solvabil¬ 
ity of regulator equations plays a key role in the study 
of nonlinear GOR. Therefore, we give the following as¬ 
sumption for the solution of regulator equations. 

Assumption 3 For i = 1 , .. .,n, there exists a smooth 
function Zi{v, fj.) with Zi(0, fj.) = 0 such that, 

= Mzi{v, p),r{v, fi)). ( 6 ) 


Under Assumption 3, letting Ui{v,iT) = — 

gi{zi(v, p),r(v, p)) and performing a coordinate trans¬ 
formation: Zi = Zi — Zi{v,p), Ci = yi — r(v,p) gives 

{zi = f^{z^,e^,v,w,p) 

\ei = gi{zi, Ci, V, w, p) -Cui-C di 

where 

fi{zi,ei,v,w,p) = fi{zi,ei,v,p) - ^^i^f2lq{v)w 
g^{zi,ei,v,w,p) = gi{zi,ei,v,p) - u,{v,p) - q{v)w 

f,{zi,ei,v,p) = Mzi,yi,p) - fi{z,{v,p),r(v,p),p) 
g^{z„ei,v,p) = gi{zi,y^,p) - g,{zi{v, p),r(v, p), p). 


Furthermore, fi{0,0,v,0, p) = 0, 5i(0,0, u, 0, = 0, 

gi{0, 0, V, p) = 0. For simplicity, the error system can be 
rewritten as 

fz = f(z,e,v,w,p) 

\e = g{z, e,v,w, p) -\- u -\- d 

where z = col(zi,..., z„), e = col(ei,..., e„), u = 
col(ui,..., u„), d = col((ii,..., dn) and /, g are suitably 
defined by (7). 

By (8), output regulation problem of the multi-agent 
system (1), (2) and (3) is transformed into a problem of 
finding a distributed control law in the form of (4) such 
that limt^.oo e{t) = 0 and the trajectory of the closed- 
loop system is well-defined for t > 0. For this purpose, 
we introduce the following assumption for the zero dy¬ 
namics (i.e., the ^-subsystem) of system (8), though we 
need not “stabilize” this subsystem (to make .2 vanish). 

Assumption 4 For any compact subset E C K." 
(h = -I- riu, -I- n^), and for i = l,...,n, there 

exists a smooth Lyapunov function 14j(-) satisfying 
aii(ll^ill) < Vz,{zi) < a 2 i(||zi||) for some smooth 
functions aii(-); o:2i{-) C /Coo, such that, for any 
{v{t),w{t),p) e E 

l(7)< -c(ii\\zi\\) +dii'yii{ei) -C 52i'y2i{w) (9) 

where 7ii(-)j are known smooth positive definite 

functions, ai{-) is a known class /Coo function satisfying 
limsups_,,g^(a“^(s^)/s) < oo, and Su, 821 are some un¬ 
known positive constants. 

Remark 3 Although the Zi-subsystem is related with v, 
the condition (9) is not restrictive since v(t) is bounded by 
Assumption 2 and fi{0,0,v,0, p) = 0. Similar assump¬ 
tions were common ly used in the study of nonlinear out¬ 
put regulation (e.g. Dona & Huand. \20lA : Iffw & Huand. 

Win:\Xu & Hu^XMm . 

3 Main Results 

In this section, we give a constructive design to solve our 
problem by investigating system (8). 

For the following non-smooth analysis, consider an 
equation x = f{x,t) with a discontinuous righthand 
side, where /: K™ x K —)• K™ is measurable and es¬ 
sentia lly locally bounded. By Proposition 3 in ICortes I 
(I 2 OO 8 II . it has a Filippov solution on [to, ti]- Let 
V: K.™ —>■ R be a locally Lipschitz continuous func¬ 
tion. V = Cl^^QY H[f]{z,t) represents the set-valued 
Lie derivative of V, wher e dV(z) denot es the Clarke’s 
generalized gradient of V (jCortes I . l2008ll . 

Then we start with a general discussion and give a sim¬ 
pler design for a special case. 
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3.1 General Discussion 

Here each agent only knows that the unknown leader’s 
input is bounded, but does not know the exact value of 
the input bound. 

To track the leader and meanwhile reject the (un¬ 
bounded) local disturbance di , we split the total control 
effort into two parts as Ui = uf + u^, where the term 
uf is to deal with di and to make yi follow yg. It is 
well-known that internal model methods were effective 
to reject modeled disturbances. Here we construct uf 
following the same technical line. 

Let P,{s) = s”ri -+Pz.np, 

be the minimal polynomial of the matrix Si and denote 
Ti = col(Ti 4 ,..., Ti^np .) with Tij = . Take 


= 


-1 

o 




~Pi,np^ 

p* 




\I/ = 

5 ^ I 


1 Olx(np.-l) 


tion, we obtain 


the leader’s unknown inputs. The gains ki and 6i are 
designe d and updated, inde pen dent of uJiiO) , t'(O), and 
I (refer lXu & Huari^. \201(\ and I Tan j. \20lJi for similar 
techniques). 

To select a proper positive function pi(-), we perform a 
coordinate transformation fji = rji — Ti — GiCi and the 
composite system of agent i becomes 

{ Zi = fi{zi,ei,v,w,pL) 
fii = Fifji + FiGiCi - Gigi{zi,ei,v,w,y.) 
ei = gi{zi,ei, v, w, p) - -I- < 

hi — Afci -t- 6i — |ep^|. 

Taking Zi = co\{zi,f}i), we have the following result. 

Lemma 1 Under Assumption 4, for any compact subset 
S C M", there are smooth Lyapunov functions Vz^{-) sat¬ 
isfying au{\\zi\\) < Vz,izi) < d 2 i(|| 2 i||)/or some smooth 
functions dii(-) and a2i(-) € ICoo with i = 1,.. .,n such 
that, for any {v{t),w{t), p) € T, and uji{t) € , 


ViiiW I(14)< + ^Ii7n(ei) + ^2i72i(w) + <^3* 


P = ^iTi, di = 'i'iTi. (10) 

The sy s tem (10) is called a steady-state generator in 
iHuand (j2004ll . which helps us reject the unwanted dis¬ 
turbance di. Since the pair (d>i, $i) is observable, there 
exists a constant matrix Gi such that Fi = -|- G'i'ki is 

Hurwitz. To asymptotically reject the disturbance di, let 

Ui = ii = Fipi + GiUi. (11) 


where 7ii(-) and 721 (•) are known smooth positive defi¬ 
nite functions, andSii, 62i, are unknown positive con¬ 
stants. 

P roof. F r om (7 ) and gi{0,Q,v, p) = 0, by Lemma 7.8 
in iHuand (I2004D . there are known positive smooth func¬ 
tions (/li(-), 4’2i{-) and an unknown constant q > 0 such 
that 


Inspired by t he ro bust adaptive control law used in 
l.liang fc Hilll (119^ . we propose the tracking control 
with a constant A > 0 for agent i (i = 1,..., n): 


WFiGiCi - Gigi{zi,ei,v,w,p)\\ 

< Ci{(j)u{zi)\\zi\\ -I- (i)2r{ei)\ei\ -I- ||w|| -I- 1). 


Ui — kiPii^Cyi'jCyi OiS^ni^Cyi') (12) 

hi — A/bj -t- Pil^Cyi’jCyi, Oi — ICpzl 

where Cyi = — yj)- For simplicity, we take 

hi{0) = 0i{O) = 0 Isee lPrald . l2003l for a similar setting). 

Then the control law for agent i can be written as 

Ui — iPi hiPii^Cyi’jeyi ^^Sgn)^^^) 

Vi = Fipi -\- GiUi (13) 

hi — A/c^ -t- pil^Cyi’^Cyi, Oi — |epj|- 


Based on Assum ption 4 and by the chan ging supply func¬ 
tions technique (jSontag &: Teel I . Il995ll , for any smooth 
function Ai(zi) > 0, there exists a smooth function 
Wziizi) satisfying d3i(||zd|) < Wz^{zi) < d4i(||zd|) for 
some smooth functions dsi)-)) d 4 i(-) G /Coo such that 


Wg.iZi) |(7)< - Ai{zi)\\Zi\\^ 

+ 4*73i(ei)ei^ -I- 54*74i(ia)||w|p 


(16) 


where 6^i, d^i are some unknown positive constants, and 
73i(’)i 74i(’) are known smooth positive functions. 


Remark 4 Different from most existing internal model 
(IM) design f o r mu l ti-agent output reg ulation (e.g. 
\Dona & Huand. \20lA : \Su & Huand. \20lA} . the control 
(13) contains two parts: the IM design itf for distur¬ 
bance rejection and the non-smooth design u\ to handle 


Let Vziizi) = KiWzp(zi) -h 2fifPiVi with a constant Ki > 
0 to be determined, where Pi is the positive definite 
solution of the Lyapunov equation PiFiPF^Pi = —lup. ■ 
Clearly, there exist smooth functionsdii(-), d 2 i(-) £ /Coo 
satisfying dudlzdl) < ^ii('Zi) < ct2 i{\\zt\\). Using (15)- 
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(16) for all {v{t), w{t), £ E, uji{t) G , we have 

VsAz^) 1(14) < -lKiAi(zi) - Acl\\Pi\\^^li{Zi)\\%\\^ 

-\%\? + Wh^i^^{e^)+^clm\\^Uei)]e^ 
+ W5,^iA^{w )+ 

Letting Ki > max{l,4cf||Pi|p}, Ai{zi) > 1 + (j)\^{zi), 
6u > max{Kij3i,4cf||Pi|p},(52i > max{Ki54i,4cf||Pi|p}, 
>_4cf||Pi|p and > [73i(ei) + <^ii(ei)]|eip, 
72i > ['yAi{w) + IJIIrulp yields the conclusion. I 

It is not hard to check that there is a Filippov solution 
of the closed-loop system consisting of (1), (2), and (13). 
Then we present one of our main results. 

Theorem 1 Under Assumptions 1-4, there exist smooth 
positive functions Pi{-) for i = l,...,n such that our 
problem is globally solved by the distributed control (13). 

Proof. The proof will be given in two steps. 

Step 1: We first constructively give smooth positive func¬ 
tions Pi{-) for i = 1,.. .,n and show the existence of the 
trajectory {zi,yi,ki,pi) for t > 0. Since 77i-subsystem is 
linear with a Hurwitz system matrix, we only have to 
show the boundedness of {zi,yi,ki). For this_purpose, 
we seek_to prove the boundedness of {z,e,k), where 
k = col(A:i,..., fcji), ki = ki — K and K > 0 is a constant 
to be determined. 

From Assumption 2, we can always find a compact sub¬ 
set E C K” containing {v{t),w{t), p) for f > 0. Based on 
Assumption 4, L emma 1 and the chan ging supply func¬ 
tions technique (ISontag fc Teel I . Il995[i , for any smooth 
function Ai{zi) > 0, there is a smooth function Wz^{zi) 
satisfying a 3 ^{\m\) < Wj,(fj) < d4i(||2j||) for some 
smooth functions Q!3i(-) and d4i(-) G /Coo such that 

lifg^(z2) 1(14) — A2(z2)||zj|| -f <^4i73i(^i)^^ 

+ S5ij4i{w)\\w\\^ -I- Se^ 

where Sa, S^i, S^i are some unknown positive constants, 
and 73i(-)j 74i(') are known smooth positive functions. 

Let Wo(z) = Ya=i 4^o(e) = e'^He. By Young’s 

inequality, the set-valued Lie derivative of Vb satisfies 

n n 

1^0 |(i4)< 3^ -I-^ ||Tij|p-l-Si-I- Pi (17) 

i=l i=l ^GdVo 

where Tij = gi{zi, ei,v, p) - - ^iGiCi and 5i = 

I]r=i(^2lk|P + ||ui(u, ^)in. Since qAO, 0, v p) = 0 and 
V is bounded, by Lemma 7.8 in lHuand (l2004li . there exist 


known positive smooth functions 4 > 2 ii') with an 

unknown positive constant Ci such that 

IIThIP < c,(0i,(z,)||z,||2 + h^ie^)\\e^\f)■ (18) 

Note that 

n n 

^^[^42732(^2) -f C2(/>22 (^2 )] 62 A ^ ^ hiCfii^Ci'jCi 
2=1 2=1 

with 6 i = S 4 i-f-Ci and a positive smooth function (e^) > 
max{73i(ei), (/) 2 i(ei)}. Recalling e„ = He, by similar ar¬ 
guments, we obtain 

n n 

2732(62) + Ci(( 2 i{Si)]s.i < E )e^2 ( 19 ) 

2=1 2=1 

for some known positive smooth functions <?!>2(-)- 

Construct a positive definite and radiahy unbounded 
Lyapunov function candidate as V{z,e,k) = Wo(z) -\- 
Vo{e) + Yn^iki - Clearly, 0^ > 0, and by (17)-(19), it 
follows 

n n 

V |(i4)<^-[Ai(z,) -Ci^H(i,)]||zi|p 

n 

- '^[2Kp{eyi) - 3 - Si^^{ey^)]eli S2 

where 

n n n 

=2 — ^1 + <55i74i(u;)||ui|p -f ^62 + XK^■ 

i—1 i—1 i—1 

For i = l,...,n, we take Ai, pi and K with Ai(zi) = 
-Ci(l)u{zi)-\-l, Pi{eyi) = 4>i{eyi) + l,K > 1 ((5*4-4). Then 

n n n 

V l(14)< - E 11^*11^ “ E “ E “2- 

i—1 i—1 i—1 

Similar to the analysis of Theorem 4.14 in lKhalill (l2002l) _. 
for any A > 0, the uniform boundedness of z, e and k 
with any solution of the closed-loop system is obtained, 
which implies the boundedness of Zi,yi, and ki. 

Step 2: Let us check e under the_controller (13). Take 
a positive definite function Vi(e,0) = Vb(e) 4- > 

where 9i = 9i — Q with 0 > 0 to be specified later. Then 

n n 

Vi 1(14) = 2^e„i[T22 - 6»iSgn(e*,i)] 4- 2^0*0^ 
i—1 i—1 

n 

< 2^e„i[T22 -0sgn(e„i)] 

i=l 
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where T 2 i = gi{zi,ei,v,w, — 'i’iGiCi — Denote Zi = co\{zi,fii), and the next lemma can be 

kiPi{evi)evi. proved in a similar way given in the last subsection. 


From the boundedness of z, e, w, w, p. and the continuity 
of related functions, there is a positive constant M such 
that ||T 2 i|| < M for all i. Letting 0 > M + 1 gives 

n n 

Vi 1 ( 14 ) <- 2 ^( 0 -M)|e„|<- 2 ^|e„|. ( 20 ) 

i=l i=l 


Lemma 2 Under Assumption 5, for any given posi¬ 
tive constant M and for i = 1,.. .,n, there is a positive 
constant M and a smooth Lyapunov function Vzj(-) 
satisfying aui\m\) < Vi,{zi) < d 2 i(||zi||) for some 
smooth functions aii(-), a 2 i{-) G /Coo, such that, for all 
{v{0),Wi{0),p) e w{t) € Wff and G 


Because Vi is radially unbounded, Ot is bounded. From 
the boundedness of e showed above, Ct is also bounded 
by (13) and (14), and then we can obtain the uniform 
continuity of V 2 = 2 X)r=i |e«l with respect to time t. 
Integrating (20) from 0 to t and taking t —>• 00 , we have 


Vi,{zi) |( 22 )< -||-2i|p + 7ii(ei) + 'J 2 i{w) + 6 

where 7 ii(-), 72 i(-) OL^e known smooth positive definite 
functions and 5 is a known positive constant. 


V 2 {e{t))dt < Fi(0). 

Recalling the Barbalat’s lemma (|Khalill . [2002tl , V 2 (e(/)) —>■ 
0 when t ^ 00 , and hence Ci converges to zero as / —)• 00 
for / = 1 ,..., n. Thus, the conclusion is obtained. I 



Then we show how the control (21) solves our problem. 

Theorem 2 Under Assumptions 1-3 and 5, our problem 
can be semi-globally solved by the distributed control ( 21 ). 

Proof. The proof is similar to that of Theorem 1. 


3.2 Special Case 


In some cases, we may know the domain of the initial con¬ 
dition (zi(0), yi(0)), u(0), Wi(0) and the bound of the un¬ 
known input and the uncertain parameter p. Of course, 
we can still use the proposed control law (13), but here 
we construct a simpler control law to solve the prob¬ 
lem based on the additional information. To this end, it 
is reasonable to introduce a new assumption to replace 
Assumption 4. 

Assumption 5 Given any compact subset S C K”, 
there exist smooth Lyapunov functions Vz^{-) satisfying 
aii(||zz||) < Vz,{zi) < a 2 i{\\zi\\) for some smooth func¬ 
tions OLii{-) and a 2 i(-) G /Coo (fori = l,...,n) such 
that, for any (v(t),w(t), p) G S, 

Vz, |(7)< -ttidl^ill) +7n(ez) +72i(w) 

where 7 ii(-) and 72 i(-) are known smooth positive defi¬ 
nite functions, and ai{-) is a known class /Coo function 
satisfying \iiLnsup,,^Q^{a~^{s"^)/s) < 00 . 

Then a new simple controller is proposed in this case: 

Ui = -'dfipi - pi{eyi)eyi - 7iSgn(e,,i) 

Vi = Fipi -\- GiUi. 


Performing a transformation Vi — Vi ~ '’'i ~ GiCi gives 


= fi{zi,ei,v,w,p) 

Vi = Fifi -\- FiGiCi - Gigi{zi, e*, v, w, p) 

Oi = gi{zi, ei,v, w, p) — "d’iVi — 'dtiGiCi -\- u^. 


Step 1: We first prove the boundedness of (£, e). Based 
on Assumption 5 and Lemma 2, we apply the chang¬ 
ing supply functions technique and obtain that, for any 
given M > 0 and any smooth function Ai{zi), there is 
a positive constant M and a smooth function Wzi{zi) 
satisfying a 3 i(||zi||) < Wzi{zi) < d4i(||.Zi|l) for some 
smooth functions d 4 i(-) G /Coo, such that, for any 

(u(0),a;i(0),^) G w{t) G and ?7,(0) G 

WiA(22)<-^^{u)m? ^ 23 ) 

+ Ai{Oi)oA + lii{w)\\w\\^ -f 

where 73 i(-), Ini') are known smooth positive functions 
and (5i is a known positive constant. 


Consider a positive definite and radially unbounded 
function V{z,e) = IFo(z) -I- Vb(e). For a given M > 0, 
by similar arguments in Theorem 1, there exist known 
smooth positive functions ifu, 4 > 2 i, such that 

IITsilP < (fiiizAlziW^ + <?i’2z(ei)||ei|P 

^ ^ (OA\ 

'^[AiAi) + (l)2i{ei)W <'^Movi)eF 

2=1 2=1 


where T 3 i = gi{zi, ei,v, p)-'d/iVi-'d/iGiei. By (17), (23), 
(24) and 7 ^ > 0, letting Ai{zi) = -cfiiiA) + I,Pz(e„i) = 
\{A{eyi)+A) yields V |( 22 )< -||i||^ - ||e,,|p-|-E; 3 , where 

n n 

S 3 = Y^iAiiw) +bl]\\w\f + nSi p)\f. 

2=1 2=1 


( 22 ) 
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Again by sim ilar techniques used in Theorem 4.14 in 
iKhalill (l2nol . we obtain the uniform boundedness of 
the trajectory {z, e) and a positive constant M, only 
depending on M and M, satisfying ||T4i|| < M with 
T 4 i = 9iizi, Ci, V, w, 9) - - '^iG^ei - Pi{eyi)evi. 

Step 2: Check e by considering the set-valued Lie deriva¬ 
tive of Vo(e) = He. Taking 7^ > M -|- 1 gives 


n 

Vo 1(22) = 2^ e^JT4* - 7iSgn(ei,i)] 

n n 

< -2^(7i - M)|e„i| < -2^ |e„i|. 

i=l i=l 

Thus, limf^+00 e(t) = 0, which implies the conclusion. I 


4 Simulations 


To illustrate our control design, we consider a multi¬ 
agent system with three non-identical followers in 
the form of high-o rder FitzHugh-Nagumo dynamics 
(|Xu fc Huand . 1201011 as follows: 

{ Xu = Xu - - X 2 i + X3i + di{t) biUi 

X2i — Xu -j- Cu X2i 

XSi — Xu -f C2i Xoi^ % — 1,2,3 

where Cii, C2i, and bi are positive constants. The lo¬ 
cal disturbances are generated by (3) with H4 = 1 -|- 
fiu Si = 0 , D 2 = [1 + M2,0], 52 = [0,1; 0,0], D 3 = 
[1 -|- /i3, 0], S 3 = [0,1; —1,0]. We aim to make Xu fol- 
l ow a reference y p = (1 + fiy)vi generated by a leader 
(iGrasmanl . llOSTH with the input p{t) unknown to the 
followers as follows: 

(vi = eo{-2vi - -yf - V 2 ) + pit) 

\w2 = Vi — V2 


where p{t) = ^ (t — r[A + (—1 )Lt + 5J is a trian¬ 

gle wave signal with period 2T and amplitude A, and 
[a;J is the largest integer not greater than x. Denote 
p, = co\{py, px, PIt P2, Ps) as the uncertain parameter 
vector. Figure 1 describes this multi-agent interaction 
topology with atj = 1, and Figure 2(a) depicts the ref¬ 
erence trajectory with ni(0) =0.1 and r;2(0) = 0. 


Withou t knowing th e exac t form of pit), the formula¬ 
tion in ISu fc Huand even with a nonlinear ex¬ 

osystem ( Chen fc Huand . 200511 fails to solve this prob¬ 
lem for our multi-agent system. Nevertheless, the prob¬ 
lem is solvable based on the formulation in Section 2. In 
fact, Assumption 2 is verified by a Lyapunov function 


V = VI + YX 2 + eoyf. Let zu = X 2 i - cu, Z 2 i = X 31 - 
C 2 i, Vi = Xu, and the plant is of the form (2) satisfying 



Fig. 1. Interaction graph Q in our example. 



Fig. 2. Trajectory of the leader and performance of the con¬ 
trol law. 

Assumption 3 with zii = (1 + py)v 2 , Z2z = “(1 + Pv)v 2 - 
Assumption 4 also holds with 14^(s) = Q!i(s) = 
7ii('S) = 72i(s) = By Theorem 1, the control 
with Piis) = s® -|- 1 can solve this problem. To reduce 
the unfavorable chattering in simulations, the sign func¬ 
tion in the proposed control (13) can be replaced by a 
saturation function as follows: 


(x/e, if|a;|<e; 

\sgn(x/e), if |a;| > e. 


With selected matrix pairs (Fi, Gi), Figure 2(b) shows 
the performance of the controller with A = 2, T = 4, 
e = 10“^, cii = C 2 i = bi = i and p G [-1, Ij®. 


5 Conclusions 

In this paper, a distributed output regulation problem 
was formulated for a class of uncertain heterogeneous 
nonlinear multi-agent systems to deal with local distur¬ 
bances and an unknown-input leader. Based on changing 
supply functions and adaptive techniques, distributed 
control laws incorporating local internal models were 
constructed. The semi-global and global results were ob¬ 
tained in two different cases. 
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